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ABSTRACT 

We use the UKIDSS Ultra-Deep Survey, the deepest degree-scale near-infrared survey to date, 
to investigate the clustering of star-forming and passive galaxies to z ^ 3.5. Our new measure- 
ments include the first determination of the clustering for passive galaxies at z > 2, which 
we achieve using a cross-correlation technique. We find that passive galaxies are the most 
strongly clustered, typically hosted by massive dark matter halos with Mhaio > 5 x 10^^ Mq 
irrespective of redshift or stellar mass. Our findings are consistent with models in which a crit- 
ical halo mass determines the transition from star-forming to passive galaxies. Star-forming 
galaxies show no strong correlation between stellar mass and halo mass, but passive galaxies 
show evidence for an anti-correlation; low-mass passive galaxies appear, on average, to be 
located in the most massive halos. These results can be understood if the termination of star 
formation is most efficient for galaxies of low stellar mass in very dense environments. 
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1 INTRODUCTION 

One of the earliest and most fundamental observational results in 
extragalactic astrophysics is t he bimodal form o f the galaxy pop- 
ulation. Ever since the time o f lHolmberg| ( ll958h . it has been clear 
that galaxies can be separated into two broad classes based on little 
more than their rest-frame optical colours. A lack of recent star- 
formation is responsible for the red colours of many galaxies, as 
the massive stars that emit light at shorter wavelengths are short 
lived. We are yet to fully determine, however, why star formation 
should cease in some galaxies yet continue in others. 

The vast majority of relatively non-star-forming (hereafter, 
passive) galaxies are also observed to have spheroidal morpholo- 
gies, while local star- forming galaxies predominantly show disks 
(e.g. iKennicuti 1 199^ . A natural approach is therefore to ask 
whether the same process responsible for morphological transfor- 
mation isEJso responsible for the termination of star formation. 
iToomrdd 19771) suggested that the merger of two similar-mass disky 
systems could result in a spheroid. A major merger, or some other 
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bulge-forming event, remains among the most favoure d explana- 
tions to date for the transition to a passive state (e.g. iBell et all 
l201lh . In this scenario a gas-rich merger or collapse of an un- 
stable disk results in \ery vigorous star formation and possibly a 
quasar phase. Either of these could have sufficient energy to expel 
the remaining gas fro m the galaxy, or heat it sufficiently to prevent 
further star formation dSilk & Reeslll998l: lFabianlll999h . Without a 
reservoir of cold gas the galaxy can no longer form stars and it is 
observed as a passive object. Episodic AGN feedback may then be 
requ ired to preverit the gas from re-coo ling onto massive galaxies 
(e.g. lBenson et alj2003| - lBest et alj20"06) . 

Theoretical considerations suggest an additional requirement 
for gala xies to end thei r star f ormation. Expanding on the early 
work of ' White & ReesI { 19781) . recent galaxy-formation models 
predict a declining rate of star formation in galaxies that depends on 
the virial t emperature of the d a rk matter halo in which the galaxy 
is loc ated i[Crotonetal.ll2006l : lOekel et al.ll2009l : iFeldmann et"al] 
l2010l : ICenll201ll) . This so called 'hot halo' model predicts that as a 
dark matter halo grows in mass it becomes able to sustain a reser- 
voir of gas with long (~ few Gyrs) cooling times. New gas infall 
is shock heated as it encounters this gas halo and is therefore not 
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immediately available for star formation. Henceforth, only a mod- 
est amount of energy input, possibly from radio AGN feedback, is 
required to maintain a high gas temperature. In this way, galaxies 
within high mass halos become starved of star-forming gas with the 
result that ongoing star formation rates become negligible. This the- 
ory explains the observed trends of downsizing: the most massive 
galaxies ar e the first with halos of sufficient mass to suppress star 
formation dCowie et al.lll996l;lKodama et alj|2004) . A clear conse- 
quence of the hot halo model is that there should be a characteristic 
halo-mass scale, above which the relative fraction of passive ob- 
jects (at fixed stellar mass) increases sharply. Current models sug- 
gest that this halo-mass scale is a few x lO^^Mp) and becomes 
maximally efficient at Mhaio ~ lO" Mq jCroton et al.ll2006l;ICenl 
I2OI 11) . This mass scale should not depend on the stellar mass of the 
galaxy sample, nor on the redshift of the objects. 

Estimating the masses of dark-matter halos that host passive 
galaxies is important to test these theoretical predictions. However, 
obtaining such measurements for individual systems is challenging, 
even in cases where excellent spectroscopic and high-resolution 
data are available. Direct halo-mass estimates (e.g. via gravitational 
lensing or velocity dispersions) are typically only available at low 
to intermediate redshifts, and only r eliable for higher-mass halos 
jMandelbaum. van de Ven. & Keeton 2009). Systematically study- 
ing the host halo masses of large populations at high redshift is 
currently beyond these techniques. Fortunately our knowledge of 
the background cosmology in which galaxies f orm and evolve is 
now sufficiently advanced jKomatsu et alJ[2oTlh such that we can 
describe the growth of the dark matter framework with reason- 
able confidence (e g. ISheth & Tonneiilll999l; IMo & White! I2OO2I : 
ISpringel et al.l2005h . This confidence has in turn allowed us to link 
galaxy samples to a typical dark matter halo mass via the use of 
their clustering statistics on linear scales (|Peebleslll980l) . 

In the low-redshift Universe, surveys such as the SDSS 
jYork et all 120001) and 2DFGRS ( Colless et alj 1200 ll) have en- 
abled studies of vast galaxy populations that had not been pos- 
sible previously. Clustering studies based upon these data rou- 
tinely find red or early type galaxy samples (those that are dom- 
inated by passive objects) cluster more strongly than an equiv- 
alent (by stel l ar ma s s or luminosity) star-forming sample (e.g . 
Norberg et id] I2OO2I; IZehavi et al.l 120051: IRoss & Brunnerl 120091: 
Zehavi et alj 201 ll) ? Moreover, contrary to star- forming samples. 



the clustering strength o f passive ga laxies does not appear to in - 
crease with luminosity ( iNorberget al . 2002; Zeh avi et al.ll2005h . 
Often, the stronger clustering of passive galaxies is attributed to 
a higher satellite fraction in t he red or passive sample dZehavi et al.l 
l2005l : lRoss & Bmnneill200^ . We return to this point in section |4] 
However, in order to fully explore the mechanisms that cause a 
galaxy to cease its star formation, we must study the population as 
the transition occurs, i.e. at higher redshift. 

At intermediate redshift (z < 1) a number of studies have 
sought to determine whether the difference between passive and 
star-forming clustering strengths continue s. Opticall y- selected sur- 
veys such as DEEP2 jPavis et alj|2003l) . VVDS ( iLe Fevre et al.l 
I2OO5I) and CFHTLS, can be used up to z ~ 1, but become heavily 
biased against red objects at higher redshift. Studies using these 
data find very simil ar results to those in the low-redshift Uni- 
verse (Meneux et aLj 200^ : ICoil et al.l2008l : lMcCracken et al.l2008l : 
ICoupon et al.B2011i) : red or passive samples are, on average, found 
in more massive dark matter halos than equivalent star-forming ob- 
jects. In fact this clustering strength depends not only on whether 
the sample is star-forming or p assive, but continuously as a func- 
tion of colour dCoil et al.ll2008l) . 



In order to use the large-scale structure to estimate halo masses 
at 2; > 1 we require selection at near-infrared wavelengths. A large 
contiguous field is clearly necessary to measure the bias at large 
scales (r ~ few Mpc) and find rare objects at high redshift, but 
equally important is excellent depth covering a wide range in wave- 
length. There have been a number of studies using a variety of 
data and colour-selection techniques to investi gate the clustering 
of passive and star-forming galaxies at z > 1 (iRoche et al.ll200a : 
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iKong et alj2"003 : lHartlev et alj2008l : rKim et alj201 ll) . However, it 
is only in recent years that survey data of the depth and breadth in 
wavelength required to construct purely mass or luminosity -limited 
samples have become available. Arguably the leading near-infrared 
survey for such work is the UKIDSS Ultra Deep Survey (UDS, see 
section |2j and these data have been used to produce two studies of 
particular r elevance to this work. 

IWilliams et alj l2009t) used the UDS data release 1 (DRl) 
to identify passive and star-forming galaxy samples based on the 
rest-fr ame colour selection criteria first introduced by IWuvts et al.l 
( l2007t) . Comparing the clustering strengths of the two samples over 
the redshift range 1 < 2; < 2, they found that passive galaxies 
are hosted by more massive halos. Breaking each sample into two 
luminosity-dependent sub-samples, they further found that there 
was no luminosity dependence on halo mass in the passive sam- 
ple. This was inteipreted as evidence for a common, or limiting, 
halo m ass for passive galax ies at this redshift. 

In lHartlevetal.l ( l2010h (hereafter, HIO) we performed a more 
extensive analysis using the third data release (DR3), over the full 
redshift range accessible from the data. We split star-forming and 
passive galaxy samples into subsamples by redshift and K-band lu- 
minosity and computed their cluste ring statistics. The res ults con- 
firmed and expanded upon those of Iwilliams et al.l ( |2009|) . finding 
a roughly constant halo mass for all passive subsamples, but clear 
evolution and generally lower halo masses for star-forming sub- 
samples. Due to the small sample of passive galaxies at 2; > 2 
however, we were unable to compute a typical halo mass for pas- 
sive galaxies at this redshift. However, the trends shown by the two 
populations suggested that passive and star-forming galaxies would 
be found in equal mass halos by 2 ~ 2.5. 

The latest UDS data are approximately 1 magnitude deeper in 
all bands compared to those used in HIO, enabling us to investigate 
the host halo masses of a stellar-mass-selected population of pas- 
sive galaxies at 2 > 2 for the first time. Furthermore, we can probe 
the environments of lower-mass passive objects (M« < 10^" Mq) 
to 2 ~ 1.5. Low-mass passive galaxies at these redshifts are ex- 
tremely faint but are a potentially important population for dis- 
criminating between galaxy formation models. Before the measure- 
ments presented in this work there were no estimates of the typical 
host halo mass for such low-mass galaxies to our knowledge. With 
these two vital pieces of information together with a consistent set 
of measurements in three stellar-mass bins across 0.3 < 2 < 3.36 
we re-visit whether there is a possible host halo-mass limit nec- 
essary to suppress star-formation in galaxies and produce passive 
objects. 

The structure of this paper is as follows: in section |2] we dis- 
cuss the data sets employed in this work and our sample construc- 
tion; section|3]details our methodology and principal results, which 
are then discussed in section|4] We conclude in section|5] Through- 
out this work, where appropriate, we adopt a ACDM cosmology 
with SI A/ = 0.3, D.A = 0.7, h= Ho/100 kms~^Mpc ~^ = 0.7 and 
(78 = 0.9. All magnitudes are given in the AB system toke & Guniil 
Il983l) unless otherwise stated. Finally, the definition of 'halo mass' 
varies in the literature; in this work we are referring to the parent 
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Figure 1. Probability-weighted photometric redshifts versus spectroscopic 
redshifts for the combined UDSz and archival redshift samples (2147 galax- 
ies). The dispersion {Az/{1 + z)) and outlier fraction are 0.031 and 3.3% 
respectively. 



halo in which a galaxy resides, as opposed to the mass of its imme- 
diate sub-halo. 



2 DATA SETS, DERIVED QUANTITIES AND SAMPLE 
DEFINITIONS 

Our sample is drawn from the K-band image of the UK IRT Infrared 
Deep Sky Survey (UKIDSS, I Lawrence et aHbOOTh . Ultra-Deep 
SurvejQ (UDS) data release 8 and matched multi-wavelength pho- 
tometry. The UDS is the deepest of the five UKIDSS sub-surveys, 
consi sting of four Wide-Field Camera (WFCAM, ICasali et ahl 
l2007h pointings, covering 0.77 square degrees in J, H and K-bands. 
Observations of the ongoing UKIDSS survey began in the spring of 
2005 and the UDS DR8 comprises all UDS data taken up to and in- 
cluding January 2010. Total on-source exposure times in the three 
near-infrared bands are 186.5, 100 and 207.5 hours for J, H and K 
respectively, reaching depths of J = 24.9, H = 24.2 and K = 24.6 
(AB magnitudes, estimated from the RMS between 2" apertures). 
The combined exposure time of these data is more than double that 
of the previous release, resulting in the deepest near-infrared im- 
ages over a single degree-scale field to date. In the case of our 
selection band (K-band), the data are more than 0.5 magnitudes 
deeper than in any comparable field. Details of the stacking proce- 
dure, mosaicing, catalogue extraction and depth estimation will be 
presented in Almaini et al. (in prep.). For this work we perform a 
catalogue extraction from the K-band image using SExtractor. We 
produced a merged source list based on two catalogues: the first 
was designed to efficiently de-blend sources and find compact ob- 
jects; while the second was optimised for extended but low-surface 
brightness galaxies. These are often referred to as 'hot' and 'cold' 
catalogues in the literature (e.g lRix et al ]l200i . 

In addition to the three UKIDSS bands, the field is covered 

^ http://www.nottingham.ac.uk/astronomyAJDS/ 
© 2012 RAS, MNRAS 0Q0.[Tlfl6l 



by comparable data from CFHT Megacam u-band, optical Subaru 
Suprime-cam data and Spitzer IRAC channels 1 and 2. The u'-band 
data reach u' = 26.75 (AB, 2" RMS) and are detailed in Foucaud 
et al. (in prep.). Deep optical data in the B, V, R, i' and z'-bands 
are taken from the Subaru-XMM Deep Survey (SXDS), achiev- 
ing depths of B = 27.6, V = 2 7.2, R = 27.0, i' = 27.0 and 
z' = 26.0 dFurusawa et al.ll2008l 5cr, 2"). Data at longer wave- 
lengths are crucial for accurate stellar-mass measurements of galax- 
ies at high redshift. The UDS Spitzer Legacy Program (SpUDS, 
PLDunlop) provides deep data in channels 1 and 2 of IRAC, as 
well as MIPS 24/im data, all of which are used during our analy- 
sis. SpUDS data reach 5o- depths of 24.2 and 24.0 (AB) at 3.6/im 
and 4.5/im respectively, while the public 24/im catalogue used here 
is limited to 300 /iJy (15cr). The co-incident area of these data sets 
after masking of bad regions and bright stars is 0.62 squar e degrees. 
Finally, we make use of existing X-ray dUeda et al.l2008h and radio 
dSimpson et alj|2006l) data to remove obvious AGN. 

Photometry was extracted in 3" apertures placed on each as- 
trom etrically alig ned image at the position of the K-band sources 
(see ISimpson et al. 2012, for further details). Due to correlated 
noise in the images that is not represented in the weight maps, mag- 
nitude uncertainties are underestimated by SExtractor. We account 
for this coiTelated noise by scaling the weight maps such that the 
uncertainty in source-free regions matches the RMS measured from 
apertures placed on the science image. 

Three of the bands (the CFHT u'-band and the two IRAC 
channels) required PSF corrections to their photometry in order to 
obtain correct colours. This correction was performed as follows. 
For the IRAC channels, the K-band image was used to estimate the 
flux lost from the aperture due to the broad IRAC PSF. The K-band 
image was smoothed to 2.0" FWHM (the FWHM of the IRAC 
images) and the ratio of counts within 3" apertures in the origi- 
nal un-smoothed image and the smoothed image was computed for 
each source. The IRAC fluxes were then multiplied by this ratio. 
In addition to correcting for the difference in PSF, this method au- 
tomatically accounts for contamination in the 3" IRAC fluxes by 
neighbouring sources. The caveat is, however, that the source and 
contaminating flux are assumed to have the same K— IRAC colour. 
We discuss this caveat further in section |4] The IRAC PSF correc- 
tion is typically of order 26 — 30% of the original flux. The u'-band 
fluxes were corrected in a similar way using the B-band image, 
however the B and u'-band PSFs are similar and so the typical cor- 
rection is less than 1% for this band. 

2.1 Photometric redshifts 

At the depths probed in our study it is essential to use photomet- 
ric redshifts (zphot)- In addition to extremely deep photometry, the 
UDS field was also the target of a unique spectroscopic survey: 
the UDSz. The UDSz is based on an ESO Large Programme tar- 
geting a large sample of galaxies (~ 3500) at Zphot > 1 with 
Kab < 23.0, plus a low-redshift control sample. The survey com- 
prised eight pointings of VIMOS in LR-Blue and LR-Red and 20 
F0RS2 masks with the GRS_300I grating. This survey has pro- 
duced ~ 1500 secure redshifts to date which are used along with 
archival redshifts (of ~ 4000 objects), after excluding objects with 
clear active galactic nuclei (AGN), to train our photometric red- 
shifts. The UDSz data will be described fully in a forthcoming pa- 
per (Almaini et al., in prep.) and a data release is expected to fol- 
low in late 2013. De tails of the archival redshifts can be found in 
ISimpson et alj ( 1201 2[) and references therein. 

Photometric redshift probability distributions were 
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Figure 2. Log(stellar mass) versus photometric redshift for our galaxy sam- 
ple (green) where the colour corresponds to the galaxy weight (paler colours 
indicate lower weights). Black circles indicate the 95% mass completeness 
limit at that redshift (see text), while the dashed line is a second order poly- 
nomial fit to these points. 



computed for each source u sing the EAZY package 
terammer. van Dokkum. & Coppil l2008h with the included 
apparent K-band magnitude prior, after correcting the ob- 
served fluxes for Galactic extincti on (following the maps of 
ISchlegel, Finkbeiner. & DavisI Il998h . The default set of six 
templates does not sufficiently represent all of our galaxies; in 
particular the u'-band flux of blue objects at high redshift is signifi- 
cantly overestimated. We therefore constructed a seve nth template, 
apply ing a small amount of SMC-like extinction jPrevot et al.l 
1 1984) to EAZY's bluest template. The differences between the 
observed and template fluxes were then used to iteratively correct 
the photometric zero points with the following adjustments relative 
to the i'-band required: u' — —0.15, 3.6fim = +0.04 and 
4.5/im = +0.05 (where a negative number means the observed 
magnitudes must be made brighter). These systematic corrections 
differ slightly from those in ICaputi et al.i ( 1201 lb . However, the 
addition of the UDSz spectra produces an overall spectroscopic 
sample which is more representative of the K-selected photometric 
sample used here than previously available samples in the field. 
The final photometric redshifts of the subset with spectra (UDSz 
plus archival) are compared with their spectroscopic redshifts in 
Figure [T] We have excluded bright X-ray sources and those with 
clear AGN signatures in their spectra, but a number of AGN are 
likely to remain in these samples. The dispersion of 2:phot — -^spcc 
for these redshifts after excluding outliers (Az/{1 + z) > 0.15, 
< 4% of objects) is Az/{1 + 2) ~ 0.031. Most previous photo- 
metric redshift sets within the field have been tested against the 
archival spectra only, which are heavily weighted towards z — 0. 
As we describe in the next section, the full probability distribution 
produced by EAZY, rather than a single probability-weighted 
redshift, is used throughout our analysis. 

2.2 Parent sample construction 

We now describe the preparation of the parent galaxy sample from 
which our passive and star-forming samples are drawn. Firstly, a 
number of cuts were made to ensure a clean but complete dataset. 
Objects were identified as stars and removed if they satisfied one of 
the following conditions: 



• For objects with K < 20, bright stellar objects are easily 
identified from a locus in light profiles defined by — . 
We make additional constraints that stars lie on the stellar locus in 
J — KvsK and in the stellar locus of a K-excess diagram (V~Jws 
J — K). Saturated objects on the K-band image (typically K < 14) 
are all assumed to be stars. 

• Objects with K > 20 are required to simultaneously lie in the 
stellar locus of a K-excess diagram, a J — K vs K diagram and a 
BzK diagram. 

These criteria are deliberately conservative as we do not wish to re- 
move any possible compact galaxies. There are 3806 objects identi- 
fied as stars in this manner. A second source of possible contamina- 
tion in the gal axy sample are sp urious sources produced by ampli- 
fier cross-talk jPye et al]|2006l) . For cross-talk produced by bright 
stars, visual identification is trivial and these are masked along with 
the stars that produce them. However, even faint objects produce 
cross-talk in deep stacks such as the UDS and these require much 
greater care to remove. The positions of cross-talk artifacts can be 
predicted, as they occur at regular intervals (128 native WFCAM 
pixels) both horizontally and vertically from the source that pro- 
duced them. Many cross-talk exhibit positive and negative pairs and 
these cases can be confirmed by cross-matching the predicted po- 
sitions with a 5a source list from an inverted image. Furthermore, 
as the optical Subaru images do not suffer from cross-talk contam- 
ination, the optical-near infrared colours can be used to select a list 
of likely artifacts. All candidate cross-talk sources from these two 
methods were then inspected in the K-band image to ensure that 
we were not removing any real objects. We find a total of 3419 
cross-talk artifact^ 

From the resulting galaxy sample a total K-band magnitude 
(SExtractor's MAG.AUTO) limit was imposed at K = 24.3 (AB). 
From simulation of sources we have determined that our catalogue 
is > 99% complete at this limit (Almaini et al. in prep.) and it 
ensures that slight variations in the depth across the UDS image do 
not affect our measurements. Following this, we also impose a pho- 
tometric redshift quality cut. Any galaxy in which the minimum 
during 2;phot fitting was greater than 11.35 was excluded. Those 
with high cannot be assigned a reliable probabilty of being at 
any redshift and will therefore contaminate our samples. Many of 
these objects contain AGN or are blended source detections, both of 
which alter the galaxies' photometry and prevent us from assigning 
them to an appropriate star-forming or passive galaxy sample. In 
addition we removed likely AGN by excluding objects associated 
with either a point-source X-ray detection or radio source brighter 
than lOOfim. Combined, these cuts remove 15% of our sample. 
We repeated our measurements with these criteria relaxed, exclud- 
ing only those with > 100 (< 1% of the sample) but found 
qualitatively identical results. We present only those mesurements 
from our cleaner sample in this paper. 

Our measurements are carried out in redshift slices, 500Mpc 
(physical) in depth. This depth corresponds to at least three times 
the photometric redshift dispersion quoted in section 12.11 Rather 
than assign each galaxy to a particular redshift, we make use of the 
full probability distribution and allow a galaxy to be represented 
multiple times. For each redshift slice, zl < z < z2, the inte- 
grated probability of a galaxy being within that redshift range is 
computed: psiicc ~ JJ^^^ p(z)ci2. Provided that psUcc > 0.02, an 
entry is made in the catalogue for that galaxy and redshift. This 

NB. A less conservative stellar catalogue was used to predict the posi- 
tions of cross-talk. 
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probability limit is imposed to reduce the computational cost of our 
measurements. For all statistical measures presented in this work, 
the integrated probability is used as a weight. Stellar masses and 
rest-frame colours (see below) are computed separately for each 
entry, at the minimum redshift within the interval. For example, 
a galaxy may have a broad peak in redshift probability, contained 
within two redshift slices, and a separate minor solution at some 
other redshift. This galaxy would be represented three times in our 
measurements, with three different values for stellar mass and three 
set s of rest-frame col ours. A similar strategy to ours was employed 
by IWake et al.l j201 l l. For simplicity we shall henceforth refer to 
each unique object and redshift as a "galaxy". 



2.3 Stellar masses, M, 

The stellar masses and rest-frame luminosities used in this work 
are measured using a multicolour stellar population fitting tech- 
nique. We fit the u'BVRi'z'JHK bands and IRAC Channels 
1 and 2 to a large grid of synthetic spectral energy distribu- 
tions (SEDs) cons tructed from the stellar population models of 
iBruzual & CharlotI i |2003) , assuming a Chabrier initial mass func- 
tion. The star-formation history is characterised by an exponen- 
tially declining model with various ages, metallicities and dust ex- 
tinctions. These models are parametrised by an age since the onset 
of SF, and by an e-folding time such that 

SFR{t) SFRo X e'i. (1) 

where the values of r can range between 0.01 and 13.7 Gyr, while 
the age of the onset of star formation ranges from 0.001 to 13.7 
Gyr. We exclude templates that are older then the age of the Uni- 
verse at the redshift of the galaxy being fit. The metallicity frac- 
tion is allowed to range from 0.0001 to 0.1, an d the dust content is 
parametrised, following Iciiarlot & Falll fcOOOl) , by Ti,, the effective 
V-band optical depth. We use values up to r„ = 2.5 with a constant 
inter-stellar medium fraction of 0.3. 

To fit the SEDs they are first scaled in the observed frame 
to the K-band magnitude of the galaxy. We then fit each scaled 
model template in the grid of SEDs to the measured photometry 
of the galaxy. We compute values for each template and select 
the best-fitting template, obtaining a coiTesponding stellar mass and 
rest-frame luminosities. 

From the catalogue of stellar masses a 95% mass com- 
pleteness lim i t as a function of redshift was found following 
IPozzetti et al.l ( l2010h . At each redshift, the galaxies within Az = 
0.05 of the target redshift were identified and sorted by their total 
(observed) K-band magnitude and the faintest 20% selected. Each 
of these objects then had its stellar mass scaled to the value that 
it would have if its magnitude were equal to our imposed limit 
(K = 24.3). The 95% mass completeness limit was taken as the 
95"^ percentile of the resulting scaled mass distribution. These 
redshift-dependent mass limits are shown in Figure |2] by small 
black circles. Also shown are a polynomial fit to the limits (dashed 
line, Miim ~ 8.27+0.87 2—0.07 z^) and our galaxy sample (green 
points). Darker points indicate greater weights. Galaxies which fall 
below Miim are not used in the subsequent analysis. 



2.4 Sample selection 

The division of our sample into passive and star-forming subsets 
could be performed in a number of ways. Here we are interested in 



those galaxies which have completed their major episode of star- 
formation and therefore show evolved stellar populations and \ery 
low residual star-formation. Our selection is based initially on the 
U, V and J-Bessel band rest-frame luminosities. These were used 
by Iwilliams et al. I i2009l) to separate evolved stellar populations 
fro m those with a gre ater fraction of recent star-fomation (see 
also lWuvts et al.ll200'^ . Due to the smal l samp le of high redshift 
(z > 2) passive galaxies, IWilliams et al.l j2009l) were only able to 
find a clear passive sequence to 2; ~ 2. Calibration of passive selec- 
tion criteria at higher redshift would require spectra which we do 
not have. However, o ur data are 1 magnitude deeper than those in 
Iwilliams etal] ( l2009h and the passive sequence does appear to con- 
tinue to higher redshift, as predicted by synt hetic evolution models . 
We therefore take a simple extension of the lWilliams et al.l ( |2009|) 
selection criteria to highe r redshift. For comple teness we list the 
colour criteria derived by Williams et al.l (|2009|) to select passive 
objects that we use in this work: 

U-V>0.88x V-J + 0.69 < 0.5) 

U-V>0.88x V-J + 0.59 (0.5 < z < 1.0) 

U-V>0.88x V-J + 0.49 (z > 1.0) 

with u — V > 1.3 and V — J < 1.6 in all cases. The remaining 
objects are assigned to our star-forming sample. Although these cri- 
teria efficiently select galaxies with predominantly old stellar pop- 
ulations, there is a possibility that the sample could still be con- 
taminated by dusty star-foming objects, edge-on discs or AGN. 
FurtheiTnore, the deep broadband photometry available across the 
electromagnetic spectrum provides much more information than is 
used in a simple UVJ selection. We therefore use the template fits 
obtained during stellar-mass determination to estimate the specific 
star-formation rate (SSFR) of each catalogue member and identify 
clear cases of dusty star-forming contaminants in the passive sam- 
ple. Objects with passive UVJ colours but SSFR > 10"* yr"^ 
were re-assigned to the star-forming sample. Furthermore, the 
24/^m data can also be used to identify red objects that harbour 
dust-enshrouded star-formation. The SFR can be estimated from 
the mid-infrared flux by 

SFR(M0yr-') = 7.8 x 10"^° L (24^m,LQ) 

dRuiopakam et al.ll201 ih . We took a simplistic approach to 24fim 
source identification, assigning each 24/im source to its closest K- 
band counterpart within 6" (though the vast majority were matched 
within 2"). For matching objects, the 24/im flux was converted to a 
luminosity (without K-correction), based solely on its redshift, and 
then to a SSFR as above. Any passive galaxy with SSFR24^m > 
7.43 X 10"" yr"^ (i.e. a stellar -mass doubling time less than the 
z = Hubble time for our cosmology) was also re-assigned to the 
star-forming sample. In total ~ 3% of passive objects selected via 
the UVJ criteria were moved to the star-forming sample, 90% of 
these were due to a 24/im detection. 

Given the small fraction of galaxies reassigned to the star- 
forming sample, it is unsurprising that relaxing these criteria has 
minimal impact on our results. In addition we have tested samples 
in which a stricter cut in the SSFR implied by the template fit is 
made. The results are, again, very similar to those presented here, 
but in making a stricter cut the risk of reassigning genuine pas- 
sive galaxies to the star-forming sample is increased. We therefore 
present the results for only the samples that we believe provide our 
cleanest separation between passive and star-forming galaxies. We 
note that the construction of our passive and star-forming samples 
differ from our earlier work (HIO). In HIO we contrasted the ex- 
tremes of the star-forming and passive samples. In this work our 
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Figure 3. The fraction of passive galaxies as a function of redshift for differ- 
ent stellar mass samples. Data are displayed only if they ai'e above the 95% 
completeness hmit, and the weighted number of galaxies in the full sample 
exceeds 100. The peak visible at 2 ^ 0.65 in the lower mass samples is due 
to a k nown cluster dominating the field at this redshift l lvan Breukelen et aU 
l200d) . Uncertainties are suppressed for clarity. 



aim is to achieve much more complete samples and we therefore 
place no requirement on star-forming galaxies to be vigorously star- 
forming. 

In Figure |3] we plot the fraction of galaxies that are assigned 
to our passive galaxy sample as a function of redshift, separated by 
stellar mass. These fractions are derived from summing the galaxy 
weights described above in each redshift range for the passive sam- 
ples and for the full sample. Samples in which the sum of galaxy 
weights in the combined passive plus star-forming sample is less 
than 100 are not represented. Fractions are unreliable for small 
samples and we do not compute clustering measurements for these 
small samples. Furthermore only samples that are above the 95% 
completeness limit are shown. There is a clear peak in the passive 
fraction of lower mass objects (log M, < 10.5 Mq) at the redshift 
correspond ing to a massive clu ster (or possible superstructure), at 
z ~ 0.65 jvan Breukelen et al . 2006). There is a possible second 
peak visible at z ~ 1.3 which is as yet unconfirmed. The uncer- 
tainties on the fractions have been suppressed for clarity. 



3 DARK MATTER HALO MASSES OF PASSIVE AND 
STAR-FORMING GALAXIES 

We now estimate the host dark matter halo masses of passive and 
star-forming sub-samples as a function of their stellar mass and red- 
shift, beginning with high redshift (2.00 < z < 3.36 - an interval 
of 1500 Mpc) passive galaxies. The typical mass of the parent ha- 
los that host a galaxy population can be estimated by measuring the 
linear bias of the population, with respect to the underlying dark 
matter density distribution. Such measurements can typically be 
achieved by computing the 2-point auto-correlation function (2pcf), 
1^(9), together with our knowledge of how the dark matter structure 
evolves with time (e.g. HIO and references therein). 

The 2pcf is the excess probability (over a random distribution) 
of finding a pair of galaxies at a given s eparation and can be esti- 
mated, using the lLandv & Szalavl h993h estimator, by 



Here, N{6) refers to the normalised number of pairs found at that 
angular separation and the subscripts DD, DR and RR represent 
pairs within the galaxy sample, pairs between the galaxy sample 
and a random catalogue and pairs within the random catalogue re- 
spectively. 

As described earlier, we use the probability of each galaxy 
lying within the desired redshift range as a weight. The pair counts 
are therefore the sum of the products of these weights over pairs 
of objects, rather than a simple pair count. The random catalogue 
must be sufficiently large such that statistical errors on 'w{6) are 
dominated by real galaxy pairs, and must follow the same selection 
function as the galaxy sample. In the case of an angular correlation 
function for a sample of even depth, this simply means that the 
same mask must be used to remove bad regions in both the galaxy 
and random catalogues. 

Finally, the constant, C, represents a coiTection for the inte- 
gral constraint, arising due to the limited field size. This constant 
depend s on the intrinsic clust ering strength and is computed, fol- 
lowing lRoche&EalesHl999l) , using the random-random pairs as 
follows: 



C - 



^NRR{e).w{0) 

T.Nrr{9) ■ 



(3) 



As the integral constraint depends on the intrinsic clustering, rather 
than the observed clustering, it must be fit simultaneously with the 
bias factor (below). Uncertainties on ■w{0) are estimated by a boot- 
strap analysis of 100 repetitions. 

In order to recover the host dark matter halo mass we must 
also know the angular correlation function of the dark matter dis- 
tribution over the same volume as our galaxy sample. For this pur- 
po se we use the non -linear dark matter power spectrum routine 
of ISmith et al] ( l2003h . computed at the mean redshift of the sam- 
ple. To obtain the non-linear dark matter angular correlation func- 
tion, this is then Fourier transformed and projected using the full, 
weighted redshift probability distribution, p(z), of the galaxy sam- 
ple (see section[2]and Figure|6}. Finally, we fit a single bias param- 
eter to the galaxy clustering measurements. 



■Wobs{0) = b X Wdni{9), 



(4) 



NDD{9)-2NDR{e) + NRR{e) 
Nrr(9) 



+ C. 



(2) 



minimising the x^- At large scales Wdm is well described by linear 
gravity theory and the non-linear correlation function, Wnon-iinear, 
is simply the same as the linear prediction, itJunoar- The assumption 
of linear bias is satisfied at these scales and a typical host halo mass 
for a galaxy sample can be estimated. 

At smaller scales non-linear effects become important, but by 
fitting the galaxy correlation function to the non-linear coiTelation 

function, using all points which satisfy lUnon- linear < 2 X UlUncar, 

we can better constrain the fit. At scales smaller than this, non- 
linear effects dominate and we caimot assume that the galaxy dis- 
tribution follows the dark matter distribution. The lower cut-off 
scale depends on the mean redshift of the galaxy sample, which 
for our 2 > 2 samples described below corresponds to 15 arc- 
seconds. It can be seen in Figure |4] that the data follow the non- 
linear DM correla tion function very well. This is unsurprising as 
ISmithetaLltoa l) used the halo model (see lCoorav & Shethl20oA 
for a review of the halo model) to determine the dark matter power 
spectrum at scales where non-linear effects become important. This 
model has previously been used to de scribe galaxy clustering statis- 
tics with considerable success (e.g. IZheng. Coil. & ZehavillioOTi: 
IRoss & Brunneill2009l) . We have verified that fits using non-linear 
scales in this way do not affect our conclusions, other than to reduce 
the uncertainties. 
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At very large separations our measurements become unreli- 
able, so we only use the data at separations up to 0.4 degrees. 
In HIO we took an alternative approach, by fitting a power law 
with slope fixed at J = —0.8 to the galaxy angular correlation 
function which was then de-projected, from which we could in- 
fer a 3-dimensional correlation length. The advantage of projecting 
the dark matter correlation function, as we do here, is that we can 
avoid the possible degeneracies introduced by simultaneously fit- 
ting a power-law slope and the integral constraint. 

3.1 The cross-correlation function 

A closely related measure to the 2pcf is the 2-point cross- 
correlation function (ccfn). The procedure is identical to that de- 
scribed for the 2pcf above, except that in the case of a ccfn we 
have two data sets, Dl and D2. The estimator is then modified as 
follows. 



NRR{e) 



(5) 



Under the assumptions that both data sets trace the same dark mat- 
ter distribution and are both linearly biased, we can use the bias 
measurement of a tracer population (6i) with the ccfn bias (612) to 
infer the bias of the second population: 



b2 = 



bi 



(6) 



Small sample sizes are one of the major difficulties that must be 
faced when using clustering to derive halo-mass estimates for pas- 
sive galaxies at high redshift. Small samples lead to large statisti- 
cal errors, dominating over other limiting factors such as sample 
variance, and consequently give very poor constraints on the halo 
mass. Due to the depth of our data we have a substantial population 
of galaxies at high redshift, though they are mostly star-forming 
objects. In using the cross-correlation statistics, we can cross a pas- 
sive population for which we want to know the bias with a full 
mass-limited tracer population and achieve much smaller statistical 
uncertainties on their bias, 62. 

Our weighted pair counts were made using a tree code 
adapt ed from Joshua Bar nes' publicly available N-body gravity 
code ( lBames&Hulll986h . Using a tree rather than a brute force 
algorithm speeds up the computation markedly, but at the cost of 
precision. We ran a number of tests and determined the optimal 
opening angle of G = 0.07. Using this value the measurements are 
almost unaffected, however the bootstrap analysis is made signifi- 
cantly faster. 




0.0 



0,001 



Figure 4. Angular two point auto and cross-correlation functions for passive 
galaxies at 2.0 < 2 < 3.36. The cross-coirelation function is computed 
using the full 95% mass-complete sample over the same redshift range. 
The cross-correlation function has been multiplied by (fe^j-f^/^tracGr) '■° 
allow direct comparison with the auto-correlation function. The two sets 
of points agree very well with each other, but using a mass-limited tracer 
population to sample the dark matter distribution provides additional galaxy 
pairs which reduces the uncertainties in the cross-correlation function. Also 
shown is the dark matter angular correlation function, multiplied by the 
passive galaxy bias squared (solid line). 
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Figure 5. Similar to Figure|4]but for star-forming galaxies at the same red- 
shift. 



3.2 Clustering of passive galaxies at 2: > 2 

In HID we suggested that if the observed clustering evolution of 
the passive galaxy population at z < 2 were to continue at higher 
redshift, then passive and star-forming galaxies would be found in 
the same mass dark matter halos at z ~ 2.5. Previously we had 
insufficient objects to obtain a meaningful measurement at this red- 
shift. However, as the latest UDS release is more than 1 magnitude 
deeper than the data used in HIO, we are now able to test this pre- 
diction. We use our passive galaxy sample, defined in Section |2] 
within the range 2.00 < z < 3.36 (a depth of 1.5 Gpc) and form a 
volume-limited sample by cutting at M* = 10'^° '* Mq. The sum of 
weights of passive objects with stellar mass, M, > 10^" * Mq, 
within this redshift range, i.e. the expected number of galaxies, 
is 1202. It is important to note that this number is dominated 



by objects with a high probability of lying in the redshift range 
2.00 < z < 3.36, i.e. with weight > 0.4 (E weight 4 = 1063). 
In other words, the clustering measurements do not simply reflect a 
larger, low-redshift population with a small probability of being at 
z>2. 

We compute measurements for both the 2pcf and ccfn, as 
detailed in the previous section. The tracer population used for 
the cross-correlation analysis was also constructed to be volume- 
limited, using all objects to the 95% completeness limit at 2 = 3.36 
within the same redshift range as the target sample (E weight = 
2450, see section l23t . These measurements are shown in Figure 
|4] where red points are the 2pcf and black points are the ccfn. The 
bias of the tracer population has been multiplied out from the ccfn 
(w' = w X bccfn/^traccr) to allow dlrcct comparisou with the 2pcf. 
Also shown is the non-linear dark matter angular correlation func- 



© 2012 RAS, MNRAS 0Q0.[Tlfl6l 



8 W. G. Hartley et al. 



passive 
star— forming 
tracer 



0.8 ■ 



0.2 - 



0.0 L 




Figure 6. Full weighted redshift probability distributions used to project the 
dark matter correlation function for our z > 2 samples. 



tion multiplied by b^, where 6 is the bias of the passive sample 
implied by the ccfn analysis (black solid line). 

As stated in section [3TT] the cross-coirelation technique relies 
on the target and tracer populations sampling the same dark mat- 
ter structure. The redshift probability distributions that are used to 
project the dark matter correlation function are shown in Figure|6] 
The three distributions are very similar, giving us confidence that 
our analysis is robust. 

Using the formalism of I Mo & Wh"it3 ( l2002h we can directly 
link the galaxy bias to a typical host dark matter halo mass. The 
bias measured for this stellar-mass-limited passive galaxy sample 
suggests they are typically hosted by halos with mass, Mhaio ~ 
1.7 X lO" M0 with bias, b = 5.23 ± 0.44E and an equivalent 
correlation length, ro = 9.9 ± 0.9 h^^Mpc. 

For comparison we plot the same measurements for star- 
forming galaxies in Figure |5] The bias measured from this sample 
is 6 = 3.76 ± 0.35, i.e. ~ 2.5<j lower than the passive galaxy sam- 
ple. This bias implies halo masses of Mhaio ^ 7.6 x lO^'^ Mq and 
correlation length, ro = 7.3 ± 0.6 h~^Mpc. 

For a stellar-mass-limited galaxy sample the host halo mass 
may also be recoved by using (sub-)halo matching, where the only 
important quatities ar e the space densities of the galaxy sample and 
their host halos (e.g. iMoster et al.ll2010l and references therein). 
We can use this technique as a consistency check of our results. 
For the combined passive and star-forming samples the number 
density of galaxies is 2.5 times the expected number density of 
halos with mass Mhaio > 7.6 x lO^'^ M©, while for the passive 
sample the number density of galaxies is ~ 5 times the expected 
number d ensity of halo s with Mhaio > 1-7 x 10^^ Mq. Using 
UDS data. lOuadrierZI|2008.) found an even greater disparity be- 
tween the number densities of distant red galaxies (DRGs) and 
their host halos at similar redshifts. However, it was later shown 




redshift 

Figure 7. Galaxy bias inferred from our cross-coiTelation analysis against 
redshift for passive and star-forming samples. In each panel the filled red 
points represent passive samples and the open blue points are star-forming 
samples. The stellar-mass range of galaxies in each panel is indicated in the 
upper left. Solid and dot-dashed lines show the bias for dark matter halos 
of various masses (in solar masses, as labelled). The dot-dashed lines are 
5 x 10^^ Mq and 5 X 10^^ Mq, upper and lower respectively. 



bv lTinker. Wechsler. & Zhend ( l2010l) that their results could be ex- 
plained by a combination of satellite galaxies and sample variance. 
We discuss the effects of multiple occupation of halos in section 
1431 



3.3 Evolution of host halo masses with redshift 

We now use the cross-correlation technique to perform a self- 
consistent analysis across all accessible redshifts (0.3 < z < 
3.36). We split our sample in redshift with bins of width 500 Mpc0 



^ The bias values quoted here and elsewhere have been corrected for a 
systematic error which is detailed in section 14.4.11 



^ This width is A z ~ 0.35 at z ~ 2 i.e. more than three times the 
expected photometric redshift dispersion 
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and in three mass intervals, lO^'^ < M/M© < 10^°", 
IqW-O ^ m/Mq < 10"^-^ and M > IO^^'^Mq for botii the 
passive and star-forming sample. As before, the tracer population 
is made up of all galaxies within the desired redshift range and with 
stellar mass greater than the 95% completeness limit at the upper 
redshift bound. Any subsample in which the lower mass limit falls 
below the 95% mass completeness limit is not considered. In Fig- 
ure |7] we display the bias determined from our cross-correlation 
method as a function of redshift for passive and star-forming ob- 
jects, with three panels representing the three bins in stellar mass. 
Red points coiTespond to our passive galaxy sample, while blue 
points represent the star-forming objects. Solid lines show how 
the clustering of dark matter halos depends on redshift at fixed 
ha lo mass (as labelled, in units of Mq) from the formalism given 
in IMo & White (2002) (and references therein). Finally, the dot- 
dashed lines are for two intermediate halo masses, 5 x 10" Mq 
and 5 x 10^^ Mq (upper and lower respectively). 

From Figure |7] a number of trends are apparent. Firstly, pas- 
sive galaxies are, on average, clearly more strongly clustered than 
their star-forming counterparts at similar stellar mass. This result 
is most apparent at lower stellar masses (M < 10"' Mq), which 
is driven in part by the structure at z ~ 0.65 mentioned earlier 
However, this behaviour remains even if we exclude the two af- 
fected redshift bins. The separation in clustering strength of pas- 
sive and star-forming objects is less clear for the more massive 
galaxies (M > W^" Mq), though still present. This difference 
is most pronounced for subsamples al z < 1 which is consistent 
with HIO, where we found a diminshing difference in clustering 
strength of passive and star-forming galaxies at higher redshift. 
Relative to HIO, however, the extra information provided by the 
cross-correlation technique allows us to greatly reduce the width 
of the redshift intervals for our subsamples. Furthermore, we are 
able to obtain halo-mass estimates for high-mass passive objects at 
z > 2 which was not previously possible. 

A second feature of these results is the mass dependence 
within samples. The least massive of the passive galaxies appear 
to be the most strongly clustered (and therefore in the most mas- 
sive dark matter halos). The bias typically decreases as the masses 
of the subsamples increase. This behaviour is c l earest at z < 1 and 
has previously been seen in iMcCracken et al.l l l2008h and HIO at 
these redshifts. Despite this mass dependence, our passive galaxy 
samples appear to have consistently high parent halo masses, with 
typical masses Mhaio ^ 5 x 10^ ^ Mq or higher We note also the 
works of IChuter et all ( l201ll) and lOuadri etaLl ( l2012b who studied 
the small-scale environments of galaxy samples in the UDS field. 
They each found that low-mass or low-luminosity passive galaxies 
are found in the densest environments at any given redshift. 

For star-forming sub-samples, galaxies of different stellar 
masses but similar redshifts do not show any clear trends. The mea- 
surements are consistent with star-forming galaxies being typically 
hosted by similar mass halos iiTespective of stellar mass. This ab- 
sence of a clear correlation for star-forming galaxies was identi- 
fied in the luminosity- selected subsamples in HIO. Similarly, the 
downsizing-like behaviour we found in HIO is also present in this 
work, with sub-samples of equal stellar mass typically hosted by 
less massive halos at low redshift . Simil ar b ehaviour has previou sly 
been identified bv lFoucaud et al.l j201(t) and lMoster et al. I j20ld) as 
well as HIO. The bias evolution of passive galaxies with redshift at 
fixed stellar mass is not as clear as that seen in the star-forming 
sample, but nevertheless appears to be present. The implications of 
these findings are considered further in the following section. 

The raw measurements presented in Figure |7] are difficult to 



fully digest. We therefore condense our results into a more intuitive 
set of plots in Figure|8] We plot the mean halo mass implied by our 
clustering results as a function of stellar mass for passive and star- 
forming subsamples in three redshift ranges: 0.5 < z < 1.0, 1.0 < 
z < 2.0 and 2.0 < z < 3.36. Within each of these intervals we do 
not ex pect strong evolution in host halo mass (see lConselice et al 
2005 and Brown et al . 2008 for studies at z < 1 and IWake et al 
201 ll for a discussion of the 1 < 2: < 2 interval) and so we can 



take a simple mean halo mass without biasing our interpretation. 
Each point in Figure|7]is converted to a typical halo mass, based on 
the measured bias, and these are then averaged over the indicated 
redshift intervals. The uncertainties on the data are the standard 
error on the mean. 

The behaviour described above can be seen more clearly in 
Figure [8] the strong clustering of low-mass passive samples, the 
lack of a stellar-mass dependence on halo mass for star-forming 
galaxies and the 'downsizing' of halo masses with redshift are all 
apparent. The averaged halo masses of all passive galaxy samples 
are greater than or equal to 5 x 10^^ M0 and within each redshift 
range the lowest-mass passive galaxies are found in the most mas- 
sive halos. This result is clearest at low redshifts, indicating a build 
up of relatively low mass passive objects in high mass dark matter 
halos. Although there is no clear stellar-mass dependence on halo 
mass for star-forming galaxies, the downsizing-like behaviour is 
clearly apparent. Typical halo masses for star-forming galaxies in 
our lowest redshift interval are ~ IO^^Mq, approaching lO" Mq 
at the highest redshifts. 



4 DISCUSSION 

4.1 A halo-mass dependent process? 

We have tested whether a halo-mass dependent effect, such as 
the formation of a hot gaseous halo, is responsible for the estab- 
lishment and build-up of the passive galaxy population. The hot 
halo favoured by theoretical models is expected to become effi- 
cient at terminating star formation when the halo mass exceeds 
Mhaio > 10^^'^Mq by shock heating newly accreted gas and 
thereby preventing the gas supply available for star formation. Pas- 
sive galaxies should therefore become much more abundant in ha- 
los of these masses. This in turn would result in a stronger cluster- 
ing signal for passive galaxies over their star-forming counterparts 
which ought to be measureable in a data set such as the UDS. 

Taken at face value, the average parent halo masses for pas- 
sive galaxy subsamples shown in Figure [8] are in agreement with a 
hot halo model and consistent with a limiting mass scale for pas- 
sivity. Passive galaxies of all masses are typically hosted by halos 
of Mhaio ^ 5 X 10^^ Mq or greater The typical host halo masses 
for star-forming galaxies show no consistent stellar-mass depen- 
dence. However, at z > 2 our star-forming samples are typically 
also found in halos with mass, Mhaio ^ 5 x 10^^ Mq. It is there- 
fore possible that the efficiency with which a hot halo can shut-off 
star-formation evolves with redshift. 

A positive stellar-mass dependence on host halo mass is ex- 
pected for galaxies at the centre of their dark matter halo, given 
that more massive dark matter halos will contain more gas avail- 
able for star-formation (at least initially). The lack of a clear mass 
dependence within the star-forming sample is likely due to a com- 
bination of multiple occupation of halos (see below) and halo-mass 
dependent transition to the passive sample, supplemented by a pos- 
sible contribution from sample variance. The inverse stellar-mass 
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Figure 8. Average halo masses implied by our bias measurements in thi'ee redshift intervals (as labelled). The data are mean halo masses for samples of fixed 
stellar mass, while the uncertainties are the standard eiTor on the mean. As in previous Figures, filled red points represent passive galaxy samples and open 
blue points are for star-forming objects. 



dependence shown by passive objects noted in the previous section 
could also result from halo occupation and a possible interpretation 
is discussed in detail below. 



4.2 Halo occupation effects 

Although our data support a halo-mass dependent process, they do 
not necessarily require it to be the hot halo that is responsible. We 
have so far neglected a number of alternative mass-dependent ef- 
fects that may be occurring. In the idealised case where each dark 
matter halo hosts only a single object, the bias allows us to recover 
a minimum halo mass. However, this is only the case for the objects 
of highest stellar mass and so for the majority of our samples we 
recover an occupation-weighted bias, 



fe= {N(M))n(M)6haio(M)/ng, 



(7) 



where n(M) is the space density of halos of mass M, (N(M)) is the 
mean number of galaxies residing in these halos, fehaio is the theo- 
retic al halo bias and n„ is the space density of the galaxy sample 
(e.g. ICoorav & Shet'i]|2002h . Satellite galaxies are typically found 
in higher-mass halos than galaxies of the same stellar mass that lie 
at the centre of their halo. The higher the fraction of satellite galax- 
ies in a sample, the further the measured typical host halo mass may 
be from the minimum mass of halos that host only central galaxies. 

We may see the influence of halo occupation in the cluster- 
ing results for the star-forming sample. Theoretical models predict 
a strong correlation between the stellar mass and host dark mat- 
ter halo mass for the central galaxy of a halo. Our measurements 
do not discriminate between central and satellite galaxies and we 
find no clear stellar-mass dependence. These observed trends can 
be explained by an increased abundance of low stellar mass satel- 
lite galaxies in high mass halos that increase the measured bias. 
However, the action of a hot halo may also help explain the lack 
of a clear trend due to the significant scatter in the stellar to halo- 
mass relation of central galaxies (e.g. 'Yang. Mo. & van den Boscl] 
[2009). halo-mass dependent quenching would preferentially re- 
move galaxies that inhabit more massive halos, leaving high stellar- 



mass galaxies that happen to be hosted by lower-mass halos in the 
star-forming sample. This process would naturally weaken any pos- 
itive stellar to halo-mass relation for star-forming galaxies. 

Given the possible influence of satellite galaxies we must con- 
sider the i mpact of satellite-related processes such as ram-pressure 
stripping jGunn & Got3l 19721: IVoUmer et al.ll200lh . Attempts have 
been made in the literature to do just this, using full halo occupation 
distribution modelling (HOD), where information on small scales 
is used t o inform (N(M)) . Most notable for this discussion is the 
work bv lTinker & Wetzell (|2010|). These authors used the angular 
clustering of IWilliams et al. I (l2009h and found that they could con- 
struct an HOD which was marginally consistent with the data {2a 
from the best fit) and in which there was no difference in minimum 
halo mass between central passive and central star-forming galax- 
ies at 1 < 2 < 2. This is the behaviour that might be expected if 
a satellite quenching process was the dominant reason for the en- 
hanced clustering of low-mass passive galaxies. As descr ibed be- 
low (section |4.3.1| ), our samples differ from those of Willia ms et al.l 
( 2009), most notably in the redshift intervals considered. It remains 
to be seen whether the minimum halo masses of star-forming and 
passive central galaxies are similar for our samples. It is possible 
that there exist passive central galaxies in low-mass halos that are 
simply much less abundant than similar-mass satellite galaxies and 
some further investigation is required. Ideally we would perform a 
consistent HOD analysis of our data in narrow redshift ranges to 
minimise redshift effects, however we defer such an analysis to a 
later paper. 

At low redshift it is possible to independently determine ap- 
proximate halo masses and identify the central galaxy of a halo. 
Central galaxies should not be subject to stripping or other satel- 
lite processes and may provide a clearer picture. More et al^ bOI ih 
used a group catalogue from the SDSS to compute average halo 
masses of passive and star-forming galaxies as a function of stellar 
mass, for central galaxies only. They find that for stellar masses, 
M* < lO^'''^ Mq, the average halo masses of passive galaxies 
are slightly greater than, but consistent with, those of the star- 
forming galaxies. However, at higher stellar masses differences 
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become apparent. The dominant central population switches from 
star-forming to passive at Mhaio ~ lO^'^ '^ Mq and there are no star- 
forming central galaxies in halos with mass, Mhaio > 10^^'* M©. 
This behaviour, and the mass scales at which it occurs, supports the 
halo-mass-scale process for ending star-formation; as does the ab- 
sence of central passive galaxies in halos of Mhaio < 10^^ Mq 
(though thi s latter point may be due also to incompleteness ef- 
fects in the iMore et al]|201ll sample). Their transition mass scale 
(Mhaio ~ 10^^'^ Mq) is lower than the scale we find for our lower 
redshift interval (0.5 < z < 1, Mhaio ~ 10^^ '^ Mq). This may 
further point towards a redshift dependence on halo-mass related 
quenching. However, given the differences in methodology, pos- 
sible satellite effects and uncertainties, further work would be re- 
quired to confirm any s uch evolution. 

Pen g et al. also used SDSS data, claiming that they 

were able to separate the environmental quenching effects from 
effects due only to stellar mass. At fixed over-density they found 
that the environmental quenching rate did not depend on the stellar 
mass of galaxies considered. Ram pressure stripping might naively 
be expected to be more efficient in low-mass galaxies as they are 
less able to gravitationally bind their gas. This finding could then 
favo ur a strangulation-like process cause d by the existence of a hot 
halo jLarson. Tinslev, & Caldweiilll98(]h . 

Recently evidence has been found for a stronger corre- 
lation between a galaxy's star-forming or passive state and 
its structural parameters (e. g. Sersic index, velocity dispersion) 
than with its stellar mass dWake. van Dokkum. & Franxl 1201 2|; 
lBelletal.|[20Tll) . This has been cited as evidence for a bulge- 
forming event being the key process in terminating star-formation, 
rather than a phenomenon due to st ellar or halo mass. However, 
IWake, Franx. & van Dokkuml ( 1201 2h find that the host halo mass 
also correlates more closely with the central galaxy's structural 
parameters (particularly velocity dispersion) than stellar mass, so 
it i s currently not clear tha t this argument must follow (though 
see lUWang, &jin3l2013l for a recen t counter-argument). Fur- 
thermore, Ivan der Wei et all ( 1201 ih and iMcLure et all (1201 2|) find 
substantial fractions of passive disks at high redshift. Nevertheless, 
the development of a bulge is likely to be an important event in a 
galaxy's evolution. Perhaps bulge-forming events result in the ex- 
pulsion of gas from a gas rich galaxy, but subsequent gas cooling 
and re-formation of a disk must prevented by a combination of en- 
vironmental effects and possibly an AGN. 



4.3 Comparison with previous work 

4.3.1 Comparison with Williams et al. (2009) 

In IWilliams et alj j2009l) the authors found constant clustering 
strength with luminosity for their passive sample. In our present 
study we find a weak anti-correlation, with stronger clustering for 
lower-mass sub-samples. In order to investigate whether this appar- 
ent discrepancy can be explained by sample select ion differences 
we ha ve tested a simple UVJ selection as used in IWilliams et al.l 
('2009). Using the same broad redshift interval (1 < z < 2) and an 
identical colour selection, but with our cross-correlation method- 
ology, template fits and stellar masses, we recover the same be- 
haviour, i.e. we find that the clustering strength does not depend 
on stellar mass within the passive sample. We suggest, therefore, 
that the build-up of passive satellites across this redshift range is 
appreciable and the use of narrow redshift intervals are important 
to identify subtle mass dependencies. 
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Figure 9. Comparison of our results with measurements from the literature. 
We have chosen two of our redshift bins that coincide with a number of 
other studies. As in previous plots, red points coiTespond to passive sub- 
samples and blue points ai'e star-forming subsamples from our work. Hori- 
zontal lines show the effective bias of dark matter halos for three different 
masses (in solar masses) as labelled. 



4.3.2 Comparison with other literature results 

The literature on clustering measurements is extensive. T h e mos t 
directly compar able works are those of 



Meneux et al 



Coil et al. (2008);'McCracken et al.'('2008');'Winiams et al.'('200' _ 
Foucaudetal. (2010); Hartley et al. (2010); Coupon et al. (201j]) 
Furusawa et al.l feouh ; IWake et al.l ( l201lh ; lLinetal.l ( I2OIII) and 
Julio et al.l ( l2012ll . several of which we have mentioned previously. 
Many of these studies quote correlation lengths, ro, while we report 
bias measurements. Following the assumptions made during com- 
putation of these ro values and our bias values, these are compara- 
ble quantities (though translating from one to the other is redshift 
de pendent). Where app ropriate we perform this conversion, using 
the lMo&Whit3 ( l2002l) models used throughout this work. 

In the following we shall consider two of our redshift in- 
tervals that coincide with the average redshifts of some recent 
results in the literature, at 2 ~ 1.3 and z ~ 2.2. Though 
we shall not cover results at z < 1 in detail, our results are 
broadly co nsistent with the measurements o f Meneux et al. (20060; 
ICoile t al. (2008); M cCracken et alj ( I2OO8I) ; iCoupon et alj ilQm 
and ljullo et aL(,2012i) . 



4.3.3 Literature comparison at z ~ 1.3 

In Figure |9] (up per panel) we plot o u r results at z ^ 1.3 a long- 
side those of | Foucaud et al.l ( I2OI0I) ; iFurusawa et al. I ( |201I[) and 
Wake et alj ( I2OII ). Each of these three studies consider only com- 
bined mass-selected samples, rather than the separate passive and 
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star-forming samples that we do here. We therefore expect the lit- 
erature measurements to lie between our points, but closer to the 
star-forming points at lower masses (as star-forming galaxies dom- 
inate by number at lower masses, see Figure|3)- Our results match 
the literature well except that for our highest-mass samples both 
passive and star-forming galaxies have slightly lower bias values 
than expected. Taking narrow redshift intervals helps to ensure that 
the underlying dark matter structure sampled by passive and star- 
forming objects is the same, but leaves our measurements slightly 
more susceptible to sample variance. It is most likely that sample 
variance is the cause of this discrepancy as it appears to affect both 
of our massive subsamples. Alternatively, it may be that there is 
a source of contamination that we have not been able to remove. 
Finally, we note that despite the fact that our redshift intervals are 
much narrower than those used for these literature results, the un- 
certainties on our bias measurements are similar or, in some cases, 
smaller due to the cross-correlation technique employed. 

4.3.4 Literature comparison at z ^ 2.2 

The lower panel of Figure |9] shows a comparison of our z ~ 2.2 
results with two set s of measureme nts from the literature. As pre- 
viously mentioned, iFurusawa et al.l ( |201 1 ) considered only mass- 
selected samples and the se points should be compared with our 
star-forming subsamples. iLin et al.l (l201lh . on the other hand , se- 
lected their objects by the BzK technique jPaddi et alj|2004h and 
then split their star-forming BzK sample by mass. Our measure- 
ments compare very well with these two literature results. We do 
not plot Lin et al.'s passive BzK sample as the mean redshift is 
lower than their other points. However, they find passive objects to 
be very strongly clustered, more so even than our results, but with 
a substantial uncertainty. 

4.4 Discussion of potential biases 

The construction of our samples can be performed in a number of 
ways and although we believe we have used the optimal selection 
method for this study, there is no single correct choice. Whether 
we use a maximum likelihood redshift, minimum or the prob- 
abilistic approach taken here, the choice could potentially affect 
the results and interpretation. Similarly, our methodology could, 
in principle, impact upon our measurements. In this sub-section we 
briefly discuss the influence of our sample construction choices and 
implicit assumptions on our results. 

4. 4. 1 Redshift distributions 

The redshift probability distribution functions, p(z), that are pro- 
duced by EAZY may not necessarily reflect the true probability 
distribution of any given galaxy. When these p{z) are weighted and 
summed, the resulting redshift distribution of a sample that is used 
to project the dark matter may therefore differ from the intrinsic 
redshift distribution. This in turn may lead to systematic errors in 
the recovered bias values. With our spectroscopic sample we are 
able to test whether this is the case and apply a correction. 

We split our spectroscopic sample into the same redshift in- 
tervals used in Figure [8] using photometric redshifts. The p{z) are 
summed in the same way as our science samples to form the pho- 
tometric n(z). This n(z) and the true spectroscopic n[z) for that 
subsample are each used to project an appropriate dark matter cor- 
relation function. The ratio of the two angular correlation functions 



that are produced in this way tells us whether we are over or un- 
derestimating the biases at these redshifts. We find that we require 
corrections to our final bias values by factors of 1.13, 1.2 and 1.4 
for z < 1,1 < z < 2 and z > 2 respectively. These are already in- 
cluded in all of our results and in the values given in the table at the 
end of this paper. In principle the required correction could depend 
not only on redshift, but also on SED and signal to noise (or equiva- 
lently, at a given redshift, stellar mass). However, our spectroscopic 
sample does not allow us to explore this possibility, and there is no 
existing sample which could be used to test such systematics for 
red galaxies at high redshift. 

4.4.2 Redshift correlations and contamination 

The second point to note is that, due to our use of a probabilistic 
redshift approach, the w(&) measurements at different redshifts are 
correlated. This is true of any binning method as photometric red- 
shift errors will scatter galaxies between redshift bins. However, 
there may be a concern that strongly clustered lower-redshift popu- 
lations are driving high bias values at high redshift. We verify that 
the vast majority of our subsamples are dominated by objects that 
have a substantial probability {p > 0.4) of being in the appropri- 
ate redshift range. The contamination by large numbers of galaxies 
that are at lower-redshifts should therefore be small. Excluding the 
small number of subsamples for which this is not strictly true from 
our analysis does not alter our conclusions. In addition to the above 
check, we repeat our measurements for samples that contain only 
galaxies withp > 0.4, finding consistent behaviour. 

There is also a possibility that contamination by stars, cross- 
talk or bright AGN that have not been removed are affecting some 
measurements. We have investigated a number of different cuts to 
our catalogues: using the photometric redshift best-fit multi- 
wavelength data to remove possible AGN, different stellar selec- 
tions and insisting upon a z-band detection for star-forming galax- 
ies to remove potential cross-talk. Each set of measurements re- 
sulted in similar behaviour across the passive and star-forming sam- 
ples. We present only those results for a fairly strict set of criteria 
(see section [Z2t with the confidence that we have removed the ma- 
jority of contaminants without biasing our results. 



4.4.3 IRAC PSF correction 

Our correction of the IRAC fluxes for the broad PSF and con- 
tamination by neighbouring sources makes the implicit assumption 
that both the source and the contaminant have the same K-3.6/im 
colour. If this assumption is false, then the IRAC fluxes could be 
inaccurate for a subset of our sources. In such cases, those objects 
may have poor fits during photometric redshift determination and 
erroneous stellar masses. The K-3.6^m colour of the full galaxy 
sample evolves with redshift, but is fairly constant at jz > 1. We test 
the impact of contaminants by comparing the redshift-dependent 
colour distribution of isolated galaxies with the broader population. 
These produce similar colour distributions, indicating that contam- 
ination is not a serious concern. Furthermore, as described above, 
we have reproduced our measurements using various cuts to re- 
move contaminants, finding very similar results. We therefore ex- 
pect that the number of objects which have close neighbours with 
very different K-3.6^m colours is small. 

We further test the impact of the IRAC PSF correction on our 
redshifts by way of a monte-carlo analysis. For two of our sam- 
ples, the 2.0 < z < 3.36 passive sample used in section \T2\ and 
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the 1.68 < z < 2.0 star-forming sample, we randomly vary the 
correction to the IRAC fluxes by between zero and twice the com- 
puted correction. The passive sample at z > 2 is important as it 
is the first such sample for which a clustering measurement has 
been made. The 1.68 < z < 2.0 star-forming sample is likely to 
be among the most difficult to achieve an accurate redshift as the 
4000A break is in the Y-band where we have no data and the Ly- 
man limit is still blue-ward of our u-band. The typical uncertainties 
in the weights that are used in our cross-correlation analysis are 
~ 4.7% and ~ 9.4% for these two samples respectively. At these 
levels we do not expect any inaccuracies in the PSF correction to 
impact upon our conclusions. 

4.4.4 Stellar mass uncertainties 

During stellar mass and star formation rate estimation we only 
considered templates with exponentially declining star formation 
rates. This approach is common in the literature, but it is becom- 
ing increasingly apparent that there is a population of objects at 
high-redshift which are be tter fit by an exponen tially increasing 
star-formation history (e.g. IPapovich et al.ll20lTI) . Our main con- 
cerns are the separation of our sample into passive and star- forming 
objects and their stellar masses. Clearly, not considering increas- 
ing star-foiTnation rates will have no impact upon their sample 
classification. However, the derived stellar masses for some high- 
redshift star-forming galaxies may be affected by our assumption, 
resulting in over-estimated stellar masses for these galaxies. Our 
most massive passive and star-forming high-redshift subsamples 
already have comparable clustering strength. If a number of the 
star-forming objects are of intrinsically lower mass, then the clus- 
tering of the true high-mass star-forming subsample should be even 
higher. It would be unexpected, but not impossible, that high-mass 
star-forming galaxies are in more massive halos than their passive 
counterparts. We anticipate, however, that any effect such as this is 
minor. 

Our stellar masses have uncertainties associated with them 
arising from photometric precision, template fitting and redshift tol- 
erance within a bin. These uncertainties are mitigated somewhat by 
our probabilistic weighting approach, but we cannot rule out cross- 
contamination between our stellar-mass bins. Th e uncertainty due 
to tem plate fitting alone has been estimated by ISolzonella et al.l 
feOlOf) for 12 bands of photometry covering the u-band to near- 
infrared at better than ~ 0.2 dex, based on simulated catalogues. 
This uncertainty reflects the variation in mass-to-light ratio that a 
galaxy may have. In addition, the luminosity of a galaxy is uncer- 
tain. Based on the breadth of our redshift bins, we can compute a 
maximum uncertainty on the luminosity of an object. Converting 
this to an error on the mass we find an average ~ 0.2 dex with 
slightly greater error, up to 0.3 dex, at low redshift and slightly 
lower, 0.1 dex, at our highest redshifts. Added in quadrature we 
find a final typical uncertainty of ^ 0.3 dex on our stellar masses, 
significantly smaller than our 0.5 dex bin width but nevertheless 
a possible cause of cross contamination. Such contamination will 
tend to reduce any clustering dependence on stellar mass and there- 
fore strengthens our conclusion regarding the stellar- mass depen- 
dent clustering of passive galaxies. 

4.4.5 Sample variance 

In our error estimate we do not take sample variance into account 
as this is best accomplished by comparing our results with indepen- 
dent fields. Our binning in redshift corresponds to only 500 Mpc in 



depth and the influence of structure is clearly visible in Figure |7] 
Within a redshift interval, however, both star-forming and passive 
objects are subject to the same structure and so their bias measure- 
ments should be similarly affected. Although our uncertainties are 
under estimated, the comparison of passive and star-forming ob- 
jects should therefore be reasonable. 



5 CONCLUSIONS 

We have used a cross-correlation function analysis to estimate the 
host halo masses of star-forming and passive galaxies extracted 
from the UKIDSS Ultra Deep Survey Xo z = 3.36. Relative to an 
auto-coiTelation function, our cross-correlation analysis provides 
much tighter constraints on the typical host halo masses of rare 
galaxy samples. This is due to the use of a much more abundant 
tracer population to sample the dark matter distribution. Passive 
galaxies were identified initially from a rest-frame colour selection 
which was then refined using template fits and 24/im data. The star- 
forming sample was made up of the remainder of the parent sample. 
Our main findings are as follows: 

We presented the first estimate of the typical host dark matter 
halo masses for a stellar-mass-selected (M, > 10^"'* M©) passive 
galaxy sample &X z > 2. This epoch is important as it is here that 
the first substantial populations of passive galaxies are found. Our 
measurements imply a typical host halo mass of Mhaio > 1.5 x 
10^'^ M0, which is higher than those of equal stellar mass z > 2 
star-forming galaxies at a 2.5cr significance level. 

We also performed a detailed analysis of our star-forming 
and passive galaxies in differential subsamples, 500Mpc in depth 
and in three stellar-mass ranges: 10®'^ < M./M© < 10^'' °, 
10^" ° < M./M© < 10^°-^ and M«/ M© > 10^° ^ We con- 
firmed many of the previous findings from lHartlev et al.l ( l20I0h but 
over a wider range in redshift and for much narrower redshift inter- 
vals. The notable results from this analysis are: 

• The mass of the host dark matter halos are clearly higher for 
passive galaxies than star-forming objects. 

• The difference in host halo masses of passive and star- 
forming galaxy samples diminishes at higher stellar mass (M, > 
10^" ° Mq) and at higher redshift. 

• All passive galaxies are hosted by high-mass halos. The mea- 
surements are consistent with a mechanism for galaxy transition to 
passivity requiring host halo masses, M ~ 5 x 10^^ M©. How- 
ever, the typical host halo masses of star-forming galaxies at 2 > 2 
are also above this limit. This may indicate a more complex situa- 
tion in which the efficiency of quenching by a hot halo is redshift 
dependent. 

• Low-mass passive galaxies (M, < 10^"" M©) are hosted 
by the most massive dark matter halos. This counter-intuitive be- 
haviour is most likely due to efficient processing of low-mass satel- 
lite galaxies in high mass halos and implies that multiple mecha- 
nisms may be involved in terminating star-formation. 

• The host halo masses of star-forming galaxies, on the other 
hand, show little dependence on stellar mass. We interpret this rel- 
ative absence of a correlation as being due to a delicate balance 
between two factors: the presence of a mild positive correlation be- 
tween stellar and halo masses for central galaxies and the greater 
abundance of low-mass satellite galaxies in more massive dark mat- 
ter halos. 

The results we have presented are an important step in deter- 
mining which process (or processes) are responsible for the termi- 
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nation of star-formation in galaxies. Our results indicate that a pro- 
cess dependent on halo mass is important for galaxies of all masses, 
but that it is likely supplemented by a process that is efficient in pro- 
cessing low-mass satellite systems, such as ram-pressure stripping. 
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Table 1. Table of values used in Figure |7] Columns are: number of objects, sum of weights (i.e. the expected number of galaxies, see section [T2\ . redshift 
interval and median redshift, stell ar-mass interval, bias , la uncertainty on the bias and the corresponding rg value and uncertainty. *Note that rg values are 
computed from the halo models of lMo & White! j2002h . and where the bias corresponds to an unphysical halo mass, a value of zero is used. 
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